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Abstract 

We consider a shear flow of a scale invariant Gaussian random velocity 
field which does not depend on the coordinates in the direction of the flow. 
We investigate a heat advection coming from a Gaussian random homo- 
geneous source. We discuss a relaxation at large time of a temperature 
distribution determined by the forced advection-diffusion equation. We 
represent the temperature correlation functions by means of the Feynman- 
Kac formula. Jensen inequalities are applied for lower and upper bounds 
on the correlation functions. We show that at finite time there is no ve- 
locity dependence of long range temperature correlations (low momentum 
asymptotics) in the direction of the flow but the equilibrium heat distri- 
bution has large distance correlations (low momentum behavior) with an 
index depending on the scaling index of the random flow and of the index 
of the random forcing. If the velocity has correlations growing with the 
distance ( a turbulent fiow) then the large distance correlations depend 
in a crucial way on the scaling index of the turbulent fiow. In such a 
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case the correlations increase in the direction of the flow and decrease in 
the direction perpendicular to the flow making the stream of heat more 
coherent. 

1 Introduction 

We investigate a heat advection in a random flow which is supposed to be "tur- 
bulent" . The turbulence is a complex phenomenon diSicult to define and avoid- 
ing a description in precise mathematical terms. The complexity of turbulence 
can be related to its dependence on the length scale relevant for undergoing 
experiments. In this paper we apply only some aspects of the turbulent flow: 
randomness of the velocity field, its self-similarity and long range correlations . 
The appearance of the turbulence should have an impact on transport phenom- 
ena described by an advection-diffusion equation of a passive scalar [1] . Such an 
equation can describe a transport of heat, a mass or some impurities. We are 
interested in the equilibrium distribution of solutions of the random advection- 
diffusion equation . The equilibrium is possible only under an external forcing 
(a heat source). We are interested in the equilibrium distribution at all scales. 
Such an equilibrium will depend on the forcing. The universality is possible 
only in the inertial range [2] [3] [4] where the external forcing should not be rel- 
evant (see ref.[5] for some recent shear fiow experiments) . Although the precise 
equilibrium distribution depends on the form of the forcing the asymptotic be- 
havior of correlation functions depends solely on the asymptotic behavior of the 
random forcing. We investigate the way the long range correlations of the fluid 
velocity influence the long range correlations of the temperature. 

We assume that there is a distinguished direction of the fluid velocity V. 
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We make a decomposition X = (x, z) e with x G -R"* and z € R^~'^ ; 
V(t, x) depends only on x G R'^ and has the non- vanishing components only in 
^D-d ^jj^ snch. a case it satisfies automatically VV = O;for physical applications 
Z) = 3 and d = 2 ov d = 1). As a typical example we could consider a fluid 
flow Vz{x,y) in the direction of the 2;-axis which does not depend on z. We can 
impose such an anisotropy of the flow by an external force R which depends 
only on x and has non-zero components solely in the z direction. So, we consider 
the Navier-Stokes equation with such a random force R 

dtV + VVV - z/AV = R 

The (0, V(x)) solution of the Navier-Stokes equation is the solution of the linear 
equation ( for the z-component) 

dtV - i^AxV = R 

( together with a zero solution for the x-componcnt) . By a proper choice of the 
external force R wc can simulate a large class of x-dependent flows. 

In sees. 2-3 wc discuss the advcction-difFusion equation, the random veloc- 
ity and a random forcing. The advection-diffusion equation can be solved by 
means of the Feynman-Kac formula. The Feynman-Kac solution has already 
been discussed by other authors [6]- [7]. These authors have been interested 
in the asymptotic behavior of the advection-diffusion equation without forcing. 
Our main interest (sees. 4-5) is in the asymptotic behavior for large time and 
distances of correlation functions of the temperature field resulting from the 
advection-diffusion equation with forcing describing the heat injection. First, in 
sec. 3 we simulate forcing by a constant gradient term in the temperature. We 
obtain a simple soluble model of advection illustrating some general features. 
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In general, we can obtain some lower and upper bounds on the correlation func- 
tions by means of the Jensen inequalities (sec. 5). For the sake of simplicity we 
concentrate on the two-point correlations. In sec. 6 we show how our methods 
can be extended to multi-point correlations. We obtain asymptotic behavior of 
the Fourier transform of the correlation functions for small and large momenta. 
We compare our methods and results (in sees. 4-6 and in the Appendix B) with 
an exactly soluble model of Kraichnan [8] [4] [9] (defined by a velocity field which 
is a white noise in time). The random advection is closely connected with a 
diffusion. In fact, under some natural assumptions random advection enforces 
diffusion [10] [11] [12] and vice versa the diffusion can be expressed as a white 
noise advection [13]. However, when we choose no diffusion (zero molecular 
diffusivity) in the initial equation of advection describing the temperature evo- 
lution then we obtain a model of advection (discussed in Appendix A) as a limit 
of the solution of the random advection-diffusion equation. The limit of zero 
molecular diffusivity has been discussed earlier in refs.[14][15]. 

In the text some positive constants arise (denoted usually as K,ci, etc.) 
which are not described at each case and are not related one to another. 



2 The advection-diffusion equation 

We consider the advection in a random velocity field V ( described in the In- 
troduction) forced by a random source / 

d^er+VV9r-Y^^r = f (1) 

where /U^ is the molecular diffusivity. If the random velocity V has correlation 
functions singular at small time then eq.(l) needs a careful interpretation. If 
the singularity of the velocity's covariance is of the form 6{t — t')D{x — x') then 



there arc two standard interpretations eitlier Ito or Stratonovitch [16] [17]. The 
difFerence between them in eq.(l)is ^D{0)V-^d. Hence, choosing one of them will 
change only the diffusion constant. We choose the Stratonovitch interpretation 
throughout the paper and also in the Appendix B. 

First, let us consider V = and / = 0. Let A'' be a (deterministic) solution 
of the heat equation 

drNr - ^ANr = (2) 
We expand 6 around the solution N of the diffusion equation 

e = T + N 

(if the mean value of V is zero then T describes fluctuations of the temperature). 
Prom eq.(l) 

drTr + WTr - yAT., = F (3) 

where 

F = f- WNr 

As the simplest example of a physical relevance we consider the mean gra- 
dient [18] [19] 

N = -gX (4) 

where g is a constant vector. The mean gradient is a stationary solution of the 
heat equation between two planes kept at fixed temperatures. For such a static 
solution 

F = / + Vg (5) 

We can see that even if / = then F is non-trivial. This is a frequent realization 
of an advection in experiments [20] [21]. In such a case the source F has the 
same distribution as the velocity. A constant mean gradient is distinguishing 

5 



a direction in space. It breaks the rotational symmetry. As a model we could 
consider g = (0, 0, .g^) and V = (0, 0, V^). 

We define the spectral measure p of the temperature T which is directly 
measurable in experiments [22] 

(T,(x,z)T,(x',z')) - (T,(x,z))(r,(x',z')) 

(6) 

= / dk(ipexp(zk(x — x') + ip(z — z'))p^(k, p) 

We have 



/dx((r,(x,p)f.(x',p')) - (f.(x,p))(f.(x',p'))) 

= <5(p + p')Pr(0,p) 



and 



(f,(x,p)7;(x,p'))-(iV(x,p))(TV(x,p)) = 5(p + p') J dkp,(k,p) 

(r^(x,z)r,(x,z)) - (T^(x,z))(r,(x,z)) = J dp J o!kp,(k,p) 

When the spectral function has singularities at low momenta then the Fourier 
transform in eq.(6) may need a careful definition in the sense of generalized 
functions. Instead of the correlation functions of ^^-(x, y) we could consider the 
structure functions 

a|'")(x,z) = ((T,(0,0) - (T,(0,0)) -T,(x,z) + (T,(x,z)))2") 

For n = 1 we have 

^r^^(x, z) = 2/ rfkdp/9r(k, p)(l — exp(ikx + ipz)) 

Qt ' scales in the same way as (TT) but has better infrared behaviour. The struc- 
ture functions ^(^") are expressed by the correlation functions of the Fourier 
transforms of T^-. 



It can bo seen that the spectral measure p of the temperature T depends on 
the spectral measure of the source / and the scaling properties of the random 
velocity field. 

3 Gaussian model of the shear flow 

We decompose the fluid velocity 

V = U + v 

into the mean value U and random fluctuations v. We assume that the velocity 
V is a Gaussian Euclidean invariant random fleld with the mean zero and 
the covariance 

(?;j(s,x)?;fe(s',x')) = Gjk{s - s',x,x') (8) 

where j,k = d + 1, D. For the sake of simplicity of the arguments we shall 
sometimes separate the time-dependence choosing G of the product form TD. 
If G is a decaying function of the distance |x — x'| then a model of the vector 
fleld V can be determined by a translation invariant G, e.g., 

Gjk{s-s',x,x') = SjkT{s-s')D{x-x') = 6jkT{s-s') j dpexp{ip{x-x'))D{p) 

(9) 

where £) is a locally integrable function. 

In a description of the turbulence we consider growing long range correla- 
tions. In such a case G cannot be translation invariant . We consider a model 
with Euclidean R"^ invariant correlation functions of v(x) — v(x'). Then 

Gjkis - s',x,x') = 6,kr{s - s')(|xr^ + |xf '3 - |x - xf (10) 

This G is positive definite if T is positive definite and < /3 < 1 (the co- 
variance (10) determines Levy's model [23] of the Brownian motion depending 
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on d-paramctcrs ). When 2/3 < 2 then the vector field v(x) does not satisfy 
the Lipschitz condition. In such a case we could expect difficulties with the 
uniqueness of the flow and the uniqueness of the solution of eq.(l) at = 0. 
Fortunately, a definition of the unique solution of eq.(l) in a weak probabilistic 
sense is possible [24] [25] even without the Lipschitz condition. 

The source / is an independent Gaussian field with the covariance 

{f{s, X, z)/(s', x', z')) = M{s - s', X - x', z - z') (11) 

We take the Fourier transform of eq.(3) in the z variable. Then, this equation 
reads 

2 2 2 

drfri^, p) + (ipV(T, x) + ^ - ^Ax)f,(x, p) = Fir, X, p) (12) 

We apply the Feynman-Kac formula [17] in order to express the solution of 
eq.(12) with the initial condition Tq € L'^{dX) in the form (the uniqueness of 
the solution is discussed in [24] [25]) 

fV(x,p) =exp(-ii^)i;[exp(-ip /;V(r-s,x + /xb(s))rfs)ro(x + Mb(r),p)]+ 

/J" dt exp(- ^'p'^^-*) )j;[exp(-ip /;-* V(r - s, x + ^ih{s))ds)F{t, x + /ib(r - i), p)] 

(13) 

In eq.(13) bj {j = 1,2, ...,rf) is the Brownian motion defined as the Gaussian 
process with the covariance [17] 

E[bj{s)bk{t)] = 6jkmin{s,t) 

We are interested in the equilibrium distribution of Tr, i.e., in the limit r ^ oo. 

When T oo and To G L'^{dX) then the first term in eq.(13) is vanishing . For 
this reason we may set Tq = from the beginning. The stationary solutions 
A'^ being harmonic functions are not square integrable in R^. Admitting such 
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functions as initial conditions wc could regain the solution TV from cq.(13) (with 
F = ) . In particular, the mean gradient (4) comes from a generalized function 
To with its support concentrated at p = 0. 

Before discussing more general correlations let us consider the constant mean 
gradient (eqs.(4)-(5)) with / = and F = gV. Then, from eq.(13) (with Tq = 0) 

TV(x, p) = 5{p)E[J^ dtgV{t, X + /xb(T - t))] (14) 

We shall see that some properties of the general advection (3) appear already 
at the level of the simple model (14). It follows from eq.(14) that 

(T,(X)) = 6{p)E[ f dtgU{t, X + /.b(T - t))] 
Jo 

and 

(f,(x,p)7;(x',p')) - (iV(x,p))(f,(x',p')) = 5{p)d{p') 
/; dt dt'E[gG{t - t', X - x' + tih{T -t)- Mb'(T - t'))g] 
We calculate the integral over time. First, if the covariance G is time-independent 

(a steady flow) then 

(r.(x,z)r.(x',zO) - (r,(x))(r,(x')) 

(loj 

= 4/x-4/dkexp(ik(x- x'))gG(k)g|k|-4(l - exp(-i^k2T))2 
Next, let us consider 

G{t-t',yi-x') = 5{t-t')D{x-x') (16) 

The covariance (16) does not have any physical foundations but the virtue of 
the assumption (16) is the solubility of the model (3) [8] (the Kraichnan model) 
in the sense that one can obtain a closed set of partial differential equations for 
the correlation fimctions (see the Appendix B). In our simplified version (14) 

(r.(x,z)r.(x',z')) - (r.(x))(r.(x')) ^^^^ 

= / rfkexp(ik(x - x'))g£'(k)g|k|-2(l - exp(-/i2k2r)) 
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If the V correlations are growing as in eq.(lO) then the expression (17) can 
be infrared divergent (especially at r = oo). In such a case we should rather 
consider 

((T,(0,0) - (T.(0,0)) -T,(x,z) + (T,(x,z)))2) 

(18) 

= 8/i-4/dk(l - exp(ikx))G(k)|k|-4(l - exp(-i^k2T))2 
In general, let 

G{t -t',x-x) = j dudkOiu, k) exp(iw(t - t') + ik(x - x')) (19) 

then 

(f.(x,p)f.(x',p')) - (TV(x,p))(r,(x',p')) = S{p)S{p') 

dt dt' J dligG{t - t', k)g exp(ik(x - x') - i/x^k^ {2^-1- t')) 
After the time integration 

(T,(x,z)r,(x',z')> - {TAX)){Tr{X')) = /dk(ia;exp(ik(x-x'))gG(u;,k)g 
(i^'lkl^ + a;2)-i|l - exp(-i/x2k2r - iujT)\^ 

(20) 

We assume that G is scale invariant 

G{ct, Ax) = c-"A-2^G(f, x) (21) 

(a + 7 < 1 if the time integral in eq.(14) is to be finite). This assumption has 
simple consequences for the heat transport. It may be not exact in mathematical 
models. As an example, for the shear flow solution of the Navier-Stokes equation 
discussed in the Introduction if Cji{(jj,\i) is the spectral function of the force 
distribution R then the spectral function of the stationary velocity distribution 
(obtained as a solution of the Navier-Stokes equation with the initial condition 
at Iq and then letting to ^ — oo) is 

Gil (a;, k) = C,,(a;, k) ((^k^)^ + a;^) (22) 
10 



We must choose a specific C in order to obtain a scale invariant G. 

We can see from eqs.(15)-(20) that at finite t the large distance behavior 
of the temperature correlations is the same as that of the velocity correlations 
because the behavior of p^- for small momenta does not change. However,if 
(v(x)v(O)) ~ jxp'' then at r = oo for a steady flow we obtain in eq.(18) 

((roo(x,z) - (T,o(x,z)) - roo(0,0) + (T,o(0,0)))2) ^ |x|2/'+4 

and for the Kraichnan model [8] 

((Too(x,z)-(r„,(x,z)) -r,o(0,0) + (T^(0,0)))2) ~ |x|2/5+2 

in eq.(17) . For a general time dependent G{t,x) of the form (19) we shall have 
the |x|^''~^"+^ behavior of the structure functions in eq.(20) if G scales 
as in eq.(21) ( 7 = — /?)• We can establish the behavior for large x — x' by 
means of a change of variables in the integrals (15)-(20) k = k|x — x'|~^ and 
u = u)\x — x'|~^ and an estimate of the remainder. Note that the long range 
correlations of the velocity field (7 < 0) lead to an increase of the temperature 
correlations. 

4 Gaussian white noise source 

In this section we consider F = / as a Gaussian random field independent of 
V. Estimates on the equilibrium distribution are simplified if the sources at 
different times are independent 

M(t - i', X - x', z - z') = 6{t - t')m{x - x', z - z') (23) 

We assume the form (23) of M as a technical simplification. This is a mathe- 
matical idealization still justified by an application of physical sources of heat 
(as heat injections independent at each time). 
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For a lower bound wo need an assumption that the dependence on x — x' 
is of the form of the Laplace transform (such an assumption includes the scale 
invariant distributions m which do not increase at large distances) either in the 
form 

m(x — x',z — z') = mi(x — x')mo(z — z') 

= /rfkrfpexp(ik(x — x') + ip(z — z'))?fii(k)mo(p) (24) 
= Jq°° dai'i{a) exp(— a|x — x'|^)mo(z — z') 
or in the Euclidean invariant way 



to(x — x', z — z') — dav{a) exp(— a(|x — x'p + |z — z'p)) 
= da J dpexp(ip(z — z')) exp(— a|x — x'p)i/(o, p) 

In eqs.(24)-(25) vi and f are non-negative functions. 
V in eq.(13) enters Tr in the form 

exp(iv(J)) 



(25) 



where 



with 



/■(J) = J du J dsv{s, u)J(s, u) 



J(s, u) = —6{s)pS{u — X — fih{T — s j) 
It follows that the expectation values of n products of are expressed by 

(exp(iv(J„))) = S{Jn) 

where S{3) is the characteristic function of the random field v. For a Gaussian 
random field 

S'(J) =exp(-^JGJ) (26) 
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Lot lis note that because of the translation invariance in the z variable of the 
source / we have a conservation of momenta 



The correlation functions (27) are expressed by the characteristic function (26) 
with J„ satisfying the condition (for n > 1) 



It follows that in the Gaussian case with the covariance (10) the part of G which 
is not translation invariant does not contribute to the correlation functions. 

We calculate the equal time expectation values of Tt ( eq.(13)with the zero 
initial condition) under the assumption that the random fields / and v are 
independent 



(t;(xi,pi) r^(x„,p, 



•n 



)) =<5(pi + ... + p„)H 



(27) 




(28) 



(f.(x,p)f,(x',p')) = <5(p + p') /;dtexp(-/x2p2(r - t)) 
£'[cxp(-?:p U(t - s, X + //b(.s))d.s) 
m(x - x' + /ib(r -t)- iih'{T - t),p)S{J2)] 



(29) 



where 



J2(u) = pe{s)6{u - X - b(r - s)) - p0{s)6{u - x' - b'(r - s)) 



For the Gaussian field (26) 



S'(J2) = exp ( - 5 /J"~* /J"~* dsds'pGo{s - s', ijh{s) - /ib(s'))p 

-| /o^~* /([~* dsds'pGo{s - s', /ub'(s) - /ub'(s'))p 

+ * /([^* dsds'pGois - s', X - x' + /xb(s) - /xb'(s'))p) 



(30) 



where Go is the translation invariant part of G. 
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If 

|m(x,p)| <if|mo|(p) (31) 
then from |<S'(J)| < 1 there follows the bound 

|(TV(x,p)TV(x',p'))| <i^(5(p + p')|mo|(p)M-'p-2(l-exp(-MVT)) (32) 

For a small p and a finite t the correlations (32) are bounded by r|mo|(p) 
whereas at r = oo by |mo|(p)p~^. 

Next, we apply the scale invariance of the Brownian motion 

h{at) = ^/ab(^) (33) 

in eq.(29).We write s = (r — t)a. Then, using the scaling properties (21) and 
(33) and denoting by Go the translation invariant part of G we can rewrite 
eqs.(29)-(30) in the form 

(2V(x, p)f,(x', pO) = S{p + p') /; exp(-M2p2(^ _ t)) 

E[exp{—ip Jq * U(r — s, x + /ib(s))rfs)m(x — x' + /x\/r — th{l) — [i^t — fb'(l),p 

exp ( - i(T - i)2-«-7 dada'pGoia - a', Mb(a) - Mb(a'))p 

-i(T - t)2-«-T dada'pGoia - a', /ib'(a) - lih'{a'))p 

+(t - t)2-"-7 dada'pGoia - a', (r - t)-3(x - x') + ^b(a) - /xb'(a'))p)] 

(34) 

For the Kraichnan model [8] T{s — s') = S{s — s') in eqs.(9)-(10),then a = 1 in 
eq.(22) and the formula (34) reads 



(7V(x, p)TV(x', p')) = 5{p + pO /; dteM-f^^PHr - *)) 

E[exp{—ip JJ" * U(r — s, x + fjh{s))ds)m{'x. — x' + fi^/r — th{l) — [x^r — fb'(l),p 

exp ( - (r - t)i-Tp£)(0)p + (r - tf-^ dapD{{T - t)-i{x - x') + nh{a) - fih' 

(35) 
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5 Jensen inequalities for the temperature corre- 
lations 

We are going to estimate the spectral measure (6)-(7) by an application of the 
Jensen inequality. We can obtain an upper bound on the correlation functions 
applying the Jensen inequality (/ duexpf > exp(/ dvf) if J du = 1) [26] to the 
time integral in eqs.(29)-(30) 

|(fV(x,p)7V(x',p'))| < 

2S{p + p') /J" dr Jo da da' exp(-/x2pV) 

^[|m(x-x' + MxA^b(l)-MV?b'(l),p)| 

exp - ir2-«-TpGo(o- - a', fjh{a) - nh{a'))p 

-ir2-"-')'pGo(a - a',Mb'(a) - fih'(a'))p 
+r2-""^pGo(<J - <J', r-3(x - x') + ^b(a) - i^h'{a'))p^] 

Let p{s, u; t, w) be the transition function for the Brownian motion to pass from 
u at time s to w at time t. Then, the expectation value (36) reads 

|(r.(x,p)f.(x',p'))| < 

25(p + p') J dudu'dv/dw' dr da da' cxp(— /i^p^r) 

p(0, 0; a', u)p{a', u; a, w)p{a, w; 1, z)p{0, 0; a', u')p{a', u'; a, W)p{a, w'; 1, z') 

|m(x - x' + /x-v/rz - yu-v/rz', p)| 

exp ^ - \r'^~"~'^pGo{a - a', jjm - /iu)p - \r'^~°'~'^pGo{a - a', /xw' - /iu')p 
+ir2-«-TpG'o(a - ct', r"* (x - x') + /xw - /iu')p 
^ 1^2-a-7pGQ(cr - a', r"* (x - x') + /iw' - Atu)p) 

(37) 
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Till now wc have kept the mean velocity U as an arbitrary non zero function. We 
can obtain a lower bound only if 



U = 



As claimed by some authors (see , e.g., the standard text-book [27]) the mean 
velocity does not play any essential role in turbulence. So, setting it equal to 
zero we do not lose much. Moreover, for the lower bound we must assume m 
of the form (24) (or (25)) with ?7io(p) > 0. Then, we can apply the Jensen 
inequality to the expectation value over the Brownian motion 

(7V(x,p)7V(x',p'))> 

^(P + P') Jo exp(-/i2p2r) davi{a) 

mo(p) oxp E 



_i^2-a-7 dada'pGoia - a',Mb'(a) - /ib'(a'))p 

+r^-"-^' Jo Jo dada'pGoicj - a', r-i(x - x') + /xb(a) - iJ.h'{a'))p 

-o|x - x' + M\/^b(l) - /xV^b'(l)|2 

For m of the form (25) the inequality (38) reads 



(38) 



Jo Jo dada'pGoia - a', fih{a) - fih{a'))p 



(39) 



(f,(x,p)f,(x',p')) > 
S{P + P') Jo dr exp(-/^2p2^-) j~ dat/{a, p) 
expE 

_ i^2-a-7 dada'pGoia - a',txh'{a) - tih'{a'))p 
+?-^"""^ Jo Jo dada'pGoia - a', r-^x - x') + nh{a) - Mb'(a'))p 
-a|x - x' + /iV^b(l) - /i\/fb'(l)|2 

The correlation functions (36)- (39) in general will essentially depend on the 
source distribution m. We consider m such that:i) mi is bounded from above 
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by a constant (cq.(31)) and in addition ii) mi(x) is decreasing like a power 217 
of |x|. From eq.(32) it follows that under the assumption (31) the limit r — + oo 
exists. We wish to estimate the correlation functions at r = oo under various 
conditions on mi (x) . Using the inequality (for ^ > 0) 

2 exp(-/i2pV->l(x-x', b)r2-"-^p2) < exp{-fi^p^r)+exp{-A{x-x', hy-'^^'^p^) 

and a change of variables in the r-integral in eqs.(36)-(37) r = t|p|~2-a-7 -^^e 
obtain ( when mi is a bounded function (31)) 

(foo(x,p)T^(x,p')) 

< 5{p + p')lmol(p) (ci^dpl - i)|p|-2 + C2^(i - IpDIpI" ^) 

where on the rhs of cq.(36) after an integral over r (which can be performed 
by a change of variables) we obtain a function |A(x — x', b)|~ ^-q-t (where b 
depends on a and a') whose expectation value is expressed by the rhs of eq.(37). 
This is an integrable function of u, w, z, u', w' and z'. Hence, it can be bounded 
by a constant C2. Under a stronger assumption that 

J (ixmi(x) < 00 (41) 

from eq.(37) we obtain in a similar way the bound 

/rfx(foo(x,p)foo(x',p')) 

(42) 

< 5(p + p')|mo|(p)(c3^(|p| - i)|p|-2 +C4^(i - IpDIpI"^) 

This is a bound on the spectral measure on the rhs of eq.(7). 

Wc wish to estimate the dependence of the correlation fmictions (34) on 
X — x' in a more explicit form. Note that if the velocity correlations arc defined 
by eq.(9) where £>(]<:) is an integrable function then on the basis of the Lebesgue 
lemma G is vanishing at large |x — x'|. In such a case the term depending on 
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X — x' in the exponential on the rhs of eq.(34) can be neglected. If m is in 
addition a slowly varying function of x — x' then 

(2V(x,p)2V(x',p')) ^ (f,(x,p)7;(x,p')) (43) 

There remains to discuss the turbulent flow (10). We are unable to prove precise 
upper bounds for large |x — x'| and general /?. However, if < 2/3 < 1 and d= 1 
then g{x.) = — Ix]^'' is a convex function [26] 

i?(^(x + y)) < ^5(x) + ^5(y) 

As a consequence 

exp ( - r2-"+'3p2r(a - a')|r-^(x - x') + fih{a) - ^ih' {a')\^f^^ 

< exp ( - ir2-"+/'p2r(a - cr')|2r-i(x - x')P^ - ^\2txh{a) - 2^b'(CT')P'^) 

Hence, under the assumption (31) (after the r-integration) the inequalities (36)- 

(37) at r = oo for < 2/3 < 1 read 

(roo(x,p)f^(x',p')) 

(44) 

<i^|p|-— |mo|(p)|x-x'|- — 
We expect the inequality (44) to hold true in general (under the assumption 
(31))for large |x— x'| because we obtain such a behavior of the two-point function 
if in a formal way we take the limit |x — x'| — > oo in eq.(34) neglecting terms of 
order |x — x'|~^. 

We discuss now the Jensen inequality (38) for the lower bound. It is suffi- 
cient to calculate the expectation value in the exponential (38). First, in the 
Kraichnan model (35) for the term —W in the exponential appearing in eq.(38) 
we obtain 

exp(-VF(x - x')) = exp f - r^-T / dkp^ (k) p 

(l - exp (ikr-i (x - x')) (l - exp (-M^k^)) j j 
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It is easy to see that 

exp(-VF(0)) = cxp ( - r^-'^ J dkpD (k) p 



, (46) 
(1 - /x-2k-2 (1 - cxp (-M'k2))) j > exp(-cri-Tp2) 

under the assumptions that pDp > |-D|p^, J rfk|£>|(k)0(|k| — i) < oo and 

/ dk\D\{k)k^0{- - |k|) < oo. 

In such a case we can take the limit r — > oo. In this hmit 

(fo,(x,p)f<^(x,p')) > S{p + p')moip)J^drJdau,{a) 
exp(-/x^p^r - cp^r^~'>' - 2a/x^r) (47) 
= mo(p)<5(p + p') /o°° rfrmi (/U\/2r) exp(— /i^p^r — cp^r^"''') 

The behavior of the integral (47) depends on the behavior of the source corre- 
lations TOi as a function of |x — x'|. If 

mi{iiy/2r)>K (48) 

then 

(foo(x, p)foo(x, p')) > (5(p+p')mo(p) (c5^(|p| -^)\p\~^+ce9{^-\p\)\pr^) 

(49) 

This lower bound coincides with the upper bound (40) (where a = 1). If mi 
satisfies a stronger condition ( < 1) 

TOi(|x|) > i^|x|-2f^ (50) 

( O > if it is to be of the form (24), i.e., ui{a) > Kop-'^) then 

(Too(x,p)foo(x,p'))>5(p+p')mo(p)(c7^(|p|--)|p|-'+'"+C8^(--|p|)|pr^ 

(51) 
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The inequality (51) results from the following estimate (for a + 7 < 1) 

drr~^ exp(— /i^p^r — cp^r^"""'"') 
= drr~^ exp(-/x2p^r - cp^r^"""''') + drr~^ exp(-/i^p^r - cp^r^"""''') 

> drr~^ exp(-(/x2p^ + cp^)r) + drr~^ exp(-(/i^p^ + cp^)r^~"~''') 

= |p|-'+'^ Jo"' exp(-(M2 + c)t)dt + |p|-^ t-^ exp(-(M2 + c)^-'^- 

(52) 

2 

where a(p) = |p| 2-0-7 and a = 1 in application to eq.(47). 

Next, we wish to estimate the behavior of the temperature correlations at 

large x — x' in the turbulent case (10) when 7 — — /3 < (if 7 > and mi is 
a bounded function then the temperature correlations are bounded from below 
and from above as functions of x— x' ,eq.(43)). First, we consider the Kraichnan 
model (35) {T{s — s') = 6{s — s') in eq.(lO)) with the mean velocity U = and 

D(k) ~ Ikj-'^+^T (53) 

The integral in eq.(45) is convergent for large k if 7 < and for small k if 
—7 < 1. We consider the model (10) with < (3 = —7 < 1. Let us change the 
integration variable in eq.(45) 

k= |x-x'|-iyrq (54) 
Then, after an estimate of the remainder 

exp(-M^(x - x')) > exp(-crp2|x - x'l^^^) (55) 

As a consequence 
(T,(x,p)2V(x',p')) 

> <5(p + p')mo(p)(M'p' + cp2|x - xf /5)-i (1 - exp(-T(M2p2 ^ ^^2^^ _ ^'\2f^))'^ 
Hence, for large |x — x'| we obtain 

(T^(x, p)Too(x', p')) > S{p + p')mo(p)c-ip-2|x - x']"'^ (56) 
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This lower bound for the Kraichnan model is the same as the upper bound (44) 
(here a = 1). 

Let us calculate the expectation value in the exponential of eq.(38) (denoted 
by —W) for the general G of eq.(19) 

W{x - x') = r2-«-T / dojdkpG (w, k) p 

(2 (i^i'k^ - t^y' (1 - (i;i2k2 _ (1 _ exp (-i^'k^ + zc^))) (57) 

- (i/|k|4+w2^-'exp(ikr-i (x-x')) |1 - exp (-i/x^k^ + iw) 

We estimate this integral at x = x' first. Similarly as in eq.(46) the scale 
invariance (21) leads to 

W{0) > cr2-«-Tp2 (58) 

if 

j dkduG {w, k) Q/U^k^ - iu^ 0{\k\ - 1) < 00 

and 

J J dkdujG (w, k) Q/u^lk]^ + ' - |k|) < 00 

Hence, under the assumption (50) on the basis of the inequalities (52) and (58) 
we have the lower bound ( generalizing that of eq.(51) to a 7^ 1) 

(foo(x,p)roo(x,p'))>(5(p+p')mo(p)(c0(|p|-^)|p|-2+2^+c'^(^-|p|)|pr^) 

(59) 

(at f2 = this lower bound coincides with the upper bound (40)). Next, if 
|x — x'l is large then for < —7 = /3 < 1 we obtain from eqs.(21) and (57) the 
lower bound 

exp(-M^(x - x')) > exp(-cpV2-o|x_x'|2/3) 

where the form of the rhs comes from a change of variables k = k'|x — x'|~^ 
and ui = co'\x — x'|~^ and an estimate of the remainder in eq.(57). 
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If wc restrict ourselves to G of the form (10) and 2/3 > 1 then we can derive 
a more precise lower bound for exp(— M^) with an application of the Holder 
inequality 

|x + y|2/3<22/3-l(|x|2/3 + |y|2/3) 

From eq. (38) and the Holder inequality we obtain after an elementary calculation 
of the expectation value over the Brownian paths 

exp(-iy(x - x')) > exp ( - Cr^-^-^p^ - cr^-^jx - x'|2'3p2^ (60) 

Hence, after a calculation of the expectation value in the exponential in eq.(38) 
the remaining r and the a integrals (from the representation (24)) in the corre- 
lation function (38) read 

/o°° / daiyi{a) exp{—W — a|x — x'p — 

— 5 Io° '^^ I davi{a) exp{—W — a|x — x'p) + i dr J daui{a) exp(— — 2ij^ra) 

where exp(— PF) is lower bounded by eq.(60). An easy estimate of this integral 
leads to the following inequality large x — x' 

(f'oo(x,p)T'oo(x',p')) 

> (5(p + p')mo(p)(Jri|pr — |x - x'|-'" + K2\v>\-^ |x - x'|-^"'") 

(61) 

where 

a =^4^ (62) 

For < 2/3 < 1 this lower boimd coincides with the upper bound (44) (derived 
for = 0). We expect that eq.(61) gives the asymptotic behavior of the two- 
point correlation function for any < 2/3 < 2 because such a behavior is a 
consequence of a formal exchange of the limit |x — x'| — > oo with the integral 
over t and the expectation value over the Brownian motion in eq.(34). 
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The lower bound (59) for small p is obtained by neglecting the |x — x'|- 
dependent term on the rhs of eq.(60). We can see from eq.(59) that if mo(p) — 
IpI""^ and mi(k) ~ |k|-''+20 ^^^^^ ^^le (ff) correlations behave as |p|-2-i'+2n 
for large momenta (short distances in the z direction), whereas the low mo- 

_ _ 2— 2n 

mentum behaviour (large distance )is|p| 2-a-7. These estimates show the 
effect of the random flow on the temperature correlations in the z direction. 
The effect on the temperature correlations in the x direction is described by 

the lower bound (61) and the upper bound (44). Again the decay of temper- 
ature correlations is determined by scaling indices of the velocity and source 
correlations. 



6 Higher order correlation functions 

Let us consider the multi-point correlation functions 

(TV(xi,Pi) TV(X2„,P2„)) = J2pairsIo '^^^■■■■dhn 

-^[Ilo-.fe) M^j - Xfe + A*bj(r - tj) - Aibfc(r - tk),Pj) 

exp ( - I /J""*' Jo~*' dsds'piGo{s - s',Xi -xi+ nhi{s) - nhi{s'))pi^] 

(63) 

where the sum is over all pairings in accordance with the Gaussian combina- 
torics. Prom (63) we have 

|(fV(xi,Pi) iV(x2„,P2„))| < T,pairs lo dti....dt2n 

U.U,k) ^(Pj + Pk)Sitj - tk) CXp(-i^2 p2(-^ _ 

^[U{j,k) - Xfc + Hhj{T - tj) - nhk{T - tk), Pj)\] < 00 

Hence, the equilibrium limit r ^ oo exists. 

If m is either of the form (24) or (25) then we can apply the Jensen inequality 
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to the expectation value in the form _B[exp /] > exp E[f]. We obtain an analogue 
of the lower bound (38). For the upper bound we apply the Jensen inequality 
to the time integral 

^M~Uo lo dsds' J J M8)Ji{s'){vUsMs'))) 

(64) 

< r-2 /; /; dsds' exp(-^ / / Jkm{s'){vk{sMs'))) 

where 

2ri 

J(s, u) = -6{s) ^ PkS{u - Xfc - /ibfe(r - s)) 

fc=i 

and the additional integral in eq. (64) is over the spatial variable u. 

We can repeat the basic estimates concerning the behavior for low z momenta 
and large x distances by means of the methods applied for the two-point cor- 
relations. First, by means of the Jensen inequalities we reduce the estimates of 
the expectation values to finite dimensional integrals. From the Jensen inequal- 
ities we can see that the correlation functions are bounded in r when r — > oo. 
Next, the results concerning the scaling behavior for 2n-point functions can be 
obtained by an introduction of spherical coordinates in the dt\...dtn integral in 
eq.(63). Then, the correlation hmctions scale in a simple way with respect to 
the temporal radius r. Let us explain such estimates in more detail for n = 2. 
Then, 

(f;(xi,pi) Tr(x4,P4)) = 5(pi +P3)(5(P2 +P4) 

Si) Jo dt2 exp(-M2p2^^ _ _ ^2p2(-^ _ ^2^) 

E[Th{xi - X3 Aibi(T - ti) - ^ih3{T - <i),pi) 

m(x2 - X4 + /ib2(r - ^2) - Mb4(T - ^2), P2) 

exp ( - E,=i,2 ir'' Io~'' dtdt'pjGoit - t',f,hj{t) - iihj{t'))pj 

+ T,j<k Jo""*' Io~*'' dtdt'pjGoit - t',Xj - Xfe + iJ,hj{t) - /ubfc(t'))Pfc)] +permut. 

(65) 
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where the sum is over perimitations of the numbers from 1 to 4 in aeeordance 
with the Gaussian combinatorics; in the sum in the exponential we set ti = 
and t2 = ti. Let t — ti = rcos9 , t — t2 = rsinO, t = ra cos 9 and t' = 
ra' sin 6. In such a case r scales in the exponential in the same way as in 
eqs.(37)-(38). The integral dtidt2 = drrdO adds an additional power of r. Under 
the assumption (31) the small p behavior of the correlation functions (65) at 
T = oo is determined by the integral 

|(f<^(xi,pi) foo(x4,P4))| 

~ |rno|(pi)|mo|(P2) drrE[exp{-r'^-°'-'^ J^jk Pj^oPk)] +permut. (66) 
- \'rho\{pi)\mo\{p2)E[\ TijkPj^oPkr^"^] +permut. 

For large distances , 7 = —(3 < and G of eq.(lO) we can expand the dependence 
on the Brownian motion in eq.(66) in powers of /u|xj — Xfe|~^. The leading order 
reads 

|(Too(xi,pi) foo(x4,P4))| ^^^^ 

~ |mo|(pi)|mo|(p2)| Z^jfcPjPfc|xj -Xfcp'^p^^ +permut. 

Note that the power describing the low p behavior in eq.(66) and large x be- 
havior in eq.(67) is twice as big as that for the two-point function (34) and (40) 
indicating the asymptotic scale invariance of the temperature ^^(x, p) at low 
momenta or large distances. This property can be extended to the 2n correlation 
functions where the scaling index is proportional to n as a consequence of the 
rfrr"~^ time integral in the spherical time coordinates. Such a behavior of the 
integrals suggests that if the velocities and the sources are scale invariant then 
the temperatures scale at large distances with the scale dimension determined 
by the two-point function. 
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7 Discussion 



The power-law behavior of turbulent velocity correlation functions and pas- 
sive scalar correlation functions in a homogeneous isotropic turbulent flow has 
been widely discussed in the literature since the basic papers of Kolmogorov [2] 
followed by Obukhov [28],Corrsin[29] and Batchelor[30] (concerning the scalar 
advection). The universal Kolmogorov | law for spectral velocity distribution as 
well as passive scalar distribution is derived by means of dimensional arguments 
(independent of any dynamical model). A statistical homogeneity and isotropy 
of the turbulence at a microscale in a sufficiently large space interval (called 
the " inertial range") is at the base of the Kolmogorov theory. Under these 
assumptions the velocity (or passive scalar) correlation functions are universal 
, i.e., independent of the source distribution rn. An experimental verification is 
not simple. Turbulent flows are usually non-homogeneous and non-isotropic at 
a macro scale. However, if a flow satisfying Kolmogorov assumptions is created 
then the spectral Kolmogorov law is satisfied in the inertial range [31]. Nev- 
ertheless, it is common for flows in nature that Kolmogorov assumptions are 
not satisfied (for some studies of such turbulent flows see [32] [33]). Even if the 
velocity is satisfying the Kolmogorov law the analogous Obukhov law for p may 
fail [32] [33] [34] . As the authors in [32] point out some problems with the verifi- 
cation of Kolmogorov's theory concerns a construction of a flow which would be 
homogeneous and isotropic in a sufficiently large inertial range ( usually bound- 
ary conditions or sources violate a global symmetry). They suggest a study of 
non-isotropic flows. 

An investigation of a general class of dynamical models of randomly forced 
Navier-Stokes and passive scalar equations is still beyond the reach of analytical 
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as well as numerical methods. A substantial progress has been achieved in the 
white noise randomly forced passive scalar (Kraichnan model) [8] [4] [9] . How- 
ever, the white noise distribution of velocities is quite unrealistic. Our main 
motivation in these studies was a derivation of the scaling behavior for veloc- 
ities which are not of the white noise type. A passive scalar in a shear flow 
independent of the coordinates in the direction of the flow was studied before 
in [6] [7] . However, these authors were interested in the anomalous free decay of 
solutions of the advection-diffusion equation. 

Our results predict a power-law of the passive scalar correlations in non- 
isotropic flows. The results depend on the source distribution m because the 
source / is present at any scale. We do not specify any inertial range in our 
model. In general, the correlations must depend on the source (for a discussion 
of random forcing see [35]). This can be seen from the detailed calculations in 
[36] [37] performed in the isotropic Kraichnan model (white noise in time) [8] [4] [9] . 
The two-point passive scalar correlations depend explicitly on the source and 
on the molecular diffusivity /i^. Only in a proper limit of the source covariance 
m and jj —^ the universal scaling law comes out. 

Before we summarize our results let us begin with simple models. First, 
consider a pure diffusion corresponding to V = 0. Then 

(f,(x,p)f,(x',p')) = 2(5(p + p')/o"drexp(-/xVr) 

E[m{yL - x' MV^b(l) - fiV^h'il), p)] 

= 2S{p + p')(27r)-°+'' /J" dr cM-f^-^P^r) 

J dudw exp{—^ — ^)m(x — x' -|- fl^/ru — /^v^w, p)] 

= n-H{p + p') / rfkexp(ik(x - x'))mi(k)mo(p)(p2 + k^)"! 

(l-exp(-A*2(p2 + k2)T)) 
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In the limit t oo 

(Too(x,z)Too(x',z')) = /"~^y rfkdpexp(ik(x-x'))exp(ip(z-z'))m(k,p)(p2+k2)"^ 

(68) 

Hence 

Poo(k,p) = /x-2(p2 + k2)-imi(k)mo(p) (69) 

Let us note that the behavior of the temperature correlations changes abruptly 
for large |z — z'| at r = oo in this simple model. At finite r it is the same 
as that of the source (say |z — z'|~'^+'^) whereas at r = oo it becomes |z — 
^ However, it ean be seen from eq.(68) that after the limit r ^ oo 
the limit /i ^ does not exist in the model without the advcetion. If we 
first take /x ^ then the subsequent limit r — > oo is linearly divergent in r. 
The strong /^-dependence of the asymptotic behavior means that this parameter 
sets a scale on time and space which determines different scaling behavior. In 
Appendix A we show that the limits ^ and r ^ oo can be interchanged 
in the model with a random advection . The correlation functions 5^^"^ in 
a non-isotropic Kraichnan model are discussed in Appendix B. The correlation 

functions (^(^"^(xi, pi, , X2ri, P2n) can be calculated exactly in the limit /t ^ 

(eq.(86)). They show no anomalous scaling (encountered in the isotropic model 
[8] [4]) as long as the points Xj are different. The scaling behavior can change 
after a transformation to the configuration space (the Fourier transform does 
not exist in the usual sense). 

Let us compare the two-point temperature correlation function (68) with the 
one in a random flow which is bounded in space and time ,i.e., G = (vv) ~ const. 
Under the assumption (31) we obtain 
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(^^.(x, p)foo(x', p')) Kd{p + p')mo(p) /q" exp(-/x2p2^ _ cpV^) (70) 

The integral (70) behaves as mo(p)p~^ for large p and as mo(p)|p|~^ for a 
small p in agreement with eq.(59) for O = a = 7 = 0. Our results of sees. 4 and 
5 give an extension of the simple observations on the temperature correlation 
functions derived in this section for a pure diffusion and for an advection by a 
uniformly bounded random flow. 

In our model (defined by the assumption that the velocity docs not depend 
on coordinates in the direction of the flow) the spectral distribution in the cor- 
responding momentum is proportional to the source distribution m as can be 
seen from eq.(34). We could consider a source / with the covariance m(x, z) 
which (approximately in a certain range as in refs.[36][37]) is independent of 
x. In such a case the spectral equilibrium distribution (6) for a pure diffu- 
sion poo(k, p) (67) is (5(k)mo(p)p~^ where the p~^ behavior comes from the 
molecular diflusivity. The temperature correlations remain independent of x 
and the limit /i ^ does not exist. A random advection is changing the be- 
havior of temperature correlations in x as well as in p. This change involves a 
non-perturbative mechanism which could not be seen in an expansion in V. It 
comes from an exponential of G in eq.(34). In particular, a steady flow bounded 
in x gives P(x>(k, p) = (5(k)|p|~^ for |p| <C i whereas for the random velocity 
growing in space with the index /3 (eq.(lO)) we have for a small k the behavior 
Poo(k, p) ~ |k| 2^ as follows from eq.(61). 

In experiments {D = 3) we could create an anisotropic flow with the Kol- 
mogorov index (10) /3=iin(i = 2or(i=l. In such a case we obtain definite 
predictions concerning the temperature distribution. This will be 



29 



|x — x'l 3 |z — z'l"^ where 

and D — d\s either 1 or 2 and there is a restriction a — /3 < 1 coming from the 
requirement of the integrabihty of the expression in the exponential of (34). 

In general, we can see from eqs.(40),(42),(44),(59) and (61) that the tur- 
bulent behavior 7 = — /3 < of the velocity field will (in comparison to pure 
diffusion) decrease the temperature correlations in the direction orthogonal to 
the flow and increase the correlations (at the fixed a)in the direction of the fiow. 
These effects contribute to the more coherent heat distribution in a turbulent 
stream. 

Appendix A: The limit — >^ 

If there is no diffusion [fi = 0) then our formulas in sees. 4-6 at finite r remain 
valid but need some interpretation. There is no expectation value over the 
Brownian motion. In such a case in some formulas (as in eqs.(34)-(35)) 7 = 0. 
Let us consider as an example the formula (34) at ij. = 

(T,(x,p)2V(x',p')) =(5(p + p')m(x-x',p) /J"rftexp(-ip /J""*rfsU(r - s,x)) 
exp ( - (r - tf-°' Jo dada'pGoia - a', 0)p 
+(r - t)2-« dada'pGoia - a', x - x')p) 

(71) 

For the Kraichnan model [8] (35) the formula (71) reads (with the Stratonovitch 
interpretation of the gradient term, see the discussion at the beginning of sec. 2). 
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(7V(x,p)2V(x',p')) =5(p + p')m(x-x',p) J^dtexp{-ip /J" *U(r-s,x)rfs) 
exp ( - (r - i)p£>o(0)p + (r - Op^o(x - x')p) 

(72) 

In the limit r — > oo and for U = we can calculate the integral over time in 
eq.(71) with the result 

(f„o(x, p)roo(x', p')) = CSip + p')m(x - x', p) 

( lo lo do'da'pGoia - cr', 0)p - dada'pGo{a - cr', x - x')pj 

(73) 

in agreement with the bounds (56) and (61). We can also calculate the higher 
order correlation functions. As an example, the four point function (65) reads 

(r,(xi,pi) f,(x4,P4)) 

= <5(P2 + P4)<5(pi + P3)m(xi - X3,pi)m(x2 - X4,P2) 

/; rfti /; dt2 cxp ( - E,=i,2 ir'' ir'' dtdt'pjGoit - t', 0)p,- 
+ Ej<fc fo~*' Io~*'' dtdt'pjGo{t - t', Xj - Xfe)pfc) + permut. 
We can obtain detailed estimates of the time integrals for any a. In some 
special cases the integrals can be explicitly calculated. In Appendix B we give 
the formula (eq.(88)) for the Kraichnan model {a = 1). For a steady flow 
(r(s) = 1 in eq.(10),a = 0) at r = oo the integration over tj gives 

{f^{xi,pi) f'oo(x4,P4)) 

= S{p2 + P4)'5(Pl + P;i)"'(xi - X3, pi)m(x2 - X4, P2) 
(^4pfp^L»o(xi - X3)i:'o(x2 - X4) 

-(PiP2)^(-Do(xi - X4) + £)o(x2 - X3) + i:)o(xi - X2) +-Do(x3 - X4))2^ ^ + permut. 

(74) 

where Do{xj — x^) = — — Xk\^^ in the model (10). 
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Appendix B:The Kraichnan model 

lfr{t — t') = 6{t — t') then we obtain a closed set of equations for the correlation 
functions 

5^")(xi,...,x„;pi,....,p„) = (TV(xi,pi)....iV(x„,p„)) (75) 

These equations have been derived by Kraichnan [8] for velocities depending 
on all coordinates. In our simplified model (9)-(10) they read (the odd order 
correlation functions are zero) 

9,5(2") ^ 1^2 ^i=2n ^ .^(2n) _ 1 + ^^(0)) E^:=f V^S?^^ 

+ T,<j,k> Pj^o(xj - Xfc)pfc5^2") + E<j,fc> -^(Pj + Pk)rn{^j - Xfc, Pj)5^^""^) {jk) 

(76) 

where Dq is the translation invariant part of D and S{jk) means that the 
coordinates Xj and Xk are lacking in <S. The term D{0) (adding to /x^) comes 
from the Stratonovitch interpretation of eq.(l).The solution of eq.(76) reads 

5^2") = exp(rA^)5^'") + /; dt exp((r - i)>t)7e5f ""'^ (77) 

If the operator A4 is strictly negative in the space Z/^(i?^'^") then the limit 
T — > oo exists and does not depend on the initial condition iSq^"' . 

We can express the solution of eq.(76) by means of the Feynman-Kac formula 
for the heat kernel 

(exp(r^)5f)(xi, ,X2„) = i;[exp( / dsW{h{s)))g{-Ki+ijbi{r), ...,X2„+/ib2„(r))] 

Jo 

(78) 

where 

Wis) = -1(^2 + Dom P| + E<,-,fe> P.-^o(x,- + nhjis) - Xfe - Mbfe(s))pfc 

(79) 
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We obtain an upper bound on the correlation functions (78) from the Jensen 
inequaUty as apphed to the time integral 

1 r 

(expr7W)ff)(xi, ,X2„) < - / dsE[exp{rW{h{s)))\g\{xi+iJ.hi{r), ...,X2n+fih2n{r))] 

r Jo 

(80) 

If g = exp h (or a superposition with positive coefficients of such functions as in 
eq.(25))then we have the lower bound from the Jensen inequality as applied to 
the expectation value 

(exp(rA^)exp/i)(xi, ,X2„) > expE[ / dsM^(b(s)))+/i(xi+;ubi(r), ...,X2„+/zb2„(r))] 

Jo 

(81) 

As an example, the formula for the two point function (in the limit r — > oo) 
with the velocity correlations defined by eq.(lO) reads 

5^^(xi,X2,pi,P2) = (5(pi +P2) /o°°drexp(-r^2pf) 

£'[exp(-pf J^dslxi - X2 +/ibi(s) - nh2{s)\'^^)m{xi +/ibi(r) - X2 - /ib2(r),pi)] 

(82) 

Then, the resulting correlation functions arc controlled from below and from 
above by the Jensen inequalities. For the lower bound (81) we obtain an explicit 
formula (using the representation (24) for mi) 

5S^(xi,X2,pi,P2) > (5(pi +P2)too(p) / diyi{a) dr exp(-r/i^pf ) 

^ (83) 
exp(— p^r^+^/i(r~2 |xi — X2I) — a|xi — X2p — 2/i^ra) 

where 

h{p) = Kp^(^+^) If' dxj^- dbb-^-^{i - (1 + 2„^Xb)-i eM-2(T+km^) 

(84) 

here K is a. positive constant. Prom eq.(84) it can easily be seen that for large 
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X — x' (small p in eq.(84)) the r-integrand in eq.(83) behaves as 

exp(— ii'rp^|xi — — a|xi — X2|^ — 

(as shown in another way in eq.(60);here a = l)leading as a consequence to the 
estimate (61) for the correlation functions. We can continue the Jensen inequal- 
ities for higher correlation functions as from eqs.(77) and (78) it follows that the 
correlation functions are again in the form of superpositions of exponentials. 

For lower order correlations a direct study of the differential equation (76) 
can be equally efficient. As an example, if I? = 3 and c? = 2 then the equation 
(76) at T = 00 (with the velocity covariance (10)) reads (here p = |xi — X2I) 

(m^ U,pd, - M V - P'P'W iP, P; M) = m(p, p) (85) 
where we defined 

S^'^\p,P\,P2\tJ') = S{pi +p2)T^'^\p,Pi\p) 

In contradistinction to the spherically symmetric case [8] eq.(85) is not explicitly 
soluble but its asymptotic solution (61) is easy to obtain. This asymptotic 
behavior is the same as the limit /x = of the solution (85) 

) (p, p- 0) = -p-2p-2/3rn(p, p) (86) 

In general, from eq.(76) the limit /z = can be obtained inductively 

5.^"^(xi,...,X2„;pi,...,p2n;0) = (|£>o(0)Ej=rPj - Y.<j,k>^3^o{^i - Xfc)pfe) 
E<ii>^(Pi +PO"^(xi -x;,pi)5^""^^(iZ;0) 

(87) 

The formulas for the asymptotic behavior (61) (a = 1) and (66) (a = 1,7 = 0) 
agree with the exact solution (87). For n = 1 the solution (82) takes the form 
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(86) whereas for n = 2 we have 



5^^(xi,...,X4;pi,...,p4;0) = (5-Do(0))Ej=iPj - E<j,fe> Pj-Do(xj - Xfe)pfej 



For the scale invariant random velocity field (10) Do(0) = and Do{xj — Xk) = 
— |xj — Xfep'^. It follows from cqs.(87)-(88) that the temperature correlation 
functions are scale invariant under scale transformations of the coordinates Xj 
as well as pj. When fj, = then the correlation functions S^"^ are singular at 
coinciding points (the limit /i — > has been studied earlier by other methods in 
[14] [15]). The bound (40) is valid for > 0. 
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